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Abstract 

Local M-operators for the classical sine-Gordon model in discrete space-time are constructed 
by convolution of the quantum trigonometric 4x4 _R- matrix with certain vectors in its "quantum" 
space. Components of the vectors are identified with r-functions of the model. This construction 
generalizes the known representation of M-operators in continuous time models in terms of Lax 
operators and classical r-matrix. 
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1 Introduction 



The r-matrix approach allows one to develop a unified treatment Q of non-linear integrable equations 
as hamiltonian systems having enough number of conserved quantities in involution. The role of the 
(classical) r-matrix is to provide a universal form of Poisson brackets for elements of Lax operators. 

An alternative approach 0], |^ consists in representing non- linear equations as 2D zero curvature 
conditions for a pair of matrix functions (called L and M operators) depending on a spectral parameter. 
Although this method avoids any reference to the hamiltonian aspects of the problem, the r-matrix is 
meaningful here, too. The alternative though less popular point of view on the r-matrix (which we are 
going to follow) is to treat it as a machine to produce M-operators from L-operators. Let us recall how 
it works. 

Let Ci{z) be a classical 2x2 X-operator on ID lattice with the periodic boundary condition Ci+m{z) — 
Ci{z); z is the spectral parameter. We assume the ultralocality - Poisson brackets of dynamical variables 
at different sites are equal to zero. The monodromy matrix 7i(z) is ordered product of L-operators along 
the lattice from the site / to Z + iV — 1: 

7i(z)= n ^^w- (1-1) 

i+A'-l>j>i 

Generating function of conserved quantities is obtained by taking trace of the monodromy matrix: 

r(z) = tr7I(z) (1.2) 

which does not depend on / due to the periodic boundary condition. Conserved quantities obtained by 
expanding logT(2;) in z are hamiltonians of commuting flows. 

All these flows admit a zero curvature representation. The generating function of corresponding 
M -operators is given by H , p| 



Mi{z-w)=T-^[w)\T^ r{-){Ti{w)®I) , (L3) 
L i 

where / is the unity matrix, r(z) is the (classical) 4x4 r-matrix. It acts in the tensor product of two 



2-dimensional spaces and, therefore, has natural block structure. In (1.3), tri means trace in the first 



space. Expanding the r.h.s. of eq. (1.3) in w, one gets A/-operators depending on the spectral parameter 
z. From the hamiltonian point of view, the zero curvature condition follows from r-matrix Poisson 
brackets for elements of the L-operator. In general, Mi{z;w) is a non-local quantity. 

A way to construct local A/-operators from the generating function is well known [|| , jl] , |3| . Suppose 
there exists a value zq of the spectral parameter such that the L-operator degenerates for all I, i.e. 
det£;(zo) = 0. This means that Ci{zo) is a 1-dimensional projector: 

A(.o)=^^, (1.4) 



where 



are two-component vectors and A; is a scalar normalization factor. (Of course this factor can be included 
in the vectors.) Components of the vectors as well as the A/ depend on dynamical variables. It is easy 
to see that Mi{z; zq) is a local quantity: 



(A 


r(z/zo) 




{Pi 


ai- 





Mi{z) = Mi{z- zo) = ' Vi ' \ (1.6) 



(note that the normalization factor cancels). The scalar product is taken in the first space only, so the 
result is a 2x2 matrix with the spectral parameter z. It generates inflnitezimal continuous time shifts. 
The zero curvature condition giving rise to non-linear integrable equations reads 

dtLi{z) = Mi+i{z)Ci{z) - Ci{z)Mi{z) . (1.7) 



preliminary draft 



2 



February 1, 2008 7:13 



The goal of this work is to extend eq. ( |1.6D to fully discretized analogues of classical 2D integrable 
models and to obtain a similar representation for Af -operators corresponding to discrete time flows. By 
the fully discretized models we mean a family of 2D partial difference integrable equations introduced 
and studied by Hirota Q-Q: discrete analogues of the KdV, the Toda chain, the sine-Gordon (SG) 
equations, etc. Our guiding principle is Miwa's interpretation of the discrete time flows, in which 
discretized integrable equations are treated as members of the same infinite hierarchy as their continuous 
counterparts. (This idea was further developed in as a general method to produce discrete soliton 
equations from continuous ones.) Although in this paper we deal with the discrete SG equation only, 
there is no doubt that the results can be more or less straightforwardly extended to other integrable 
models on 2D lattice. 



Let us outline the results. A formula for discrete M-operators of the type (1.6) does exist. Re- 
markably, the structure of the formula remains the same, with the classical r-matrix being substituted 
by quantum i?-matrix. Specifically, we show that the following representation of discrete M-operators 
holds: 

Mi{z) = ^ ' I 4 \(3i)=ai\Pi) (1.8) 

(here and below cti, (72, are the Pauli matrices). In the r.h.s., R{z) is a quantum 4x4 i?-matrix to be 
specified below with the "quantum" parameter q related to the discrete time lattice spacing. We stress 
that vectors |q;/) and in this formula are the same as in eq. ( |1.6| ). 



The M-operator (1.8) generates shifts of a discrete time variable to. The discrete zero curvature 
condition 

Ml+l,miz)Cl^rniz) (z)A1/^„(z) (1.9) 

(with the same L-operator as before) gives rise to non-linear integrable equations in discrete space-time 
listed above. Under the condition that these equations are satisfied a similar i?-matrix representation 
for the L-operator itself is valid: 

{(3i\R(-Hz/zo)\ai) 

Here i?*-~'(z) is another quantum i?-matrix which differs from R(^z^ by Drinfeld's twist. 

As soon as the quantum i?-matrix comes into play, we c a n use customary conventions of the algebraic 
Bethe ansatz approach ll|] to represent formulas (1.8), (1.10|) in the figure: 



(/3|i?(z)|a) 



The scalar product is taken in the " quantum" (vertical) space of the i?-matrix, so the result is a 2 x 2 
matrix in its auxiliary (horizontal) space. 

The change of dynamical variables to the pair of vectors ja;), plays a key role in the construction. 
In the papers Q on exact lattice regularization of integrable models components of these vectors were 
expresed in terms of canonical variables of the model. Those formulas looked complicated and were 
hardly considered as something instructive. Here we add a remark which may help to understand their 
meaning. Taking into account equations of motion of the completely discretized model derived from the 
discrete zero curvature condition (1.9), we show that (suitably normalized) components of the vectors 
I a/), |/?;) are t -functions. The r-function is known to be one of the most fundamental objects of the 
theory (see e.g. (l2)). 

In this paper we explain the construction for the simplest example ~ the lattice SG model. In the 
literature, there are two lattice versions of the classical SG model: Hirota's discrete analogue of the SG 
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equation on a space-time lattice ^] and the model on a space lattice with continuous time introduced 
by Izergin and Korepin Q and further studied in [l^]. Actually, they are closely connected with each 
other and, furthermore, each of the two models can be better understood with the help of the other 



one. They have common L-operator (which can be represented in the form (1.10)). The Af -operators 
are different: in the latter case the r.h.s. of eq. (1.6) with the classical trigonometric r-matrix generates 
continuous time derivatives while in the former case discrete time shifts are generated by (^^. The 
J?-matrix entering the r.h.s. of this formula is in this case the simplest trigonometric solution of the 
Yang-Baxter equation (the i?-matrix of the XXZ spin chain). At last, we demonstrate how eq. ( |1.6| ) is 
reproduced from eq. (|1.8|) in the continuous time limit. 



2 Lattice SG models with discrete and continuous time 

By the discrete SG model on the space-time lattice we mean the Faddeev-Volkov version [H of Hirota's 
discrete analogue m of the classical SG equation: 



iyip{u, V + l)ip{u + l,v + 1) — h'ip{u, v)'4'{u + 1, v) 
^■ijj{u + 1, v)%jj{u + 1, u + 1) — fJ'ipiu, v)ip{u, V + 1) 



(2.1) 



This is a non-linear partial difference equation for a function tp of the two variables where /i, v are 
constant parameters. We call u,v chiral space-time coordinates. It was shown |p^ , [l5| that eq. (2.1) 
contains both SG and KdV equations as results of different continuum limits. 

Let 



■ C={u,v+1) 



A={u,v) 



D=(«-|-l,ii+l) 



B=(u+l,v) 



be an elementary cell of the m, u-lattice; in this notation eq. (2.1) reads 



(2.2) 



The non- linear equation (2.2) can be represented |14| as a discrete zero curvature condition. Set 

(j)(u,v) ='>l)^{u,v) (2.3) 
and consider the chiral L-operator Q, |16| 



L 




(2.4) 



Equation (2.2) is equivalent to the zero curvature condition 

Ld^b{z; v)Lb^a{z; h) = Ld^c{z; ijl)Lc^a{z; v) (2.5) 

on the chiral space-time lattice. 

The chiral i-operators are building blocks for more complicated ones. For our purpose we need to 
pass from chiral to the "direct" space-time coordinates / = ^{u-\- v), m = ^{u — v). (From now on we 
refer to / (resp., m) as discrete space (resp., discrete time) coordinate.) To do that, consider "composite" 
operators which generate shifts along the diagonals A ^ D and C ^ B respectively: 



Cd^a{z) = z ^Ld^c{z; fJ-)Lc^A{z;v), 



(2.6) 
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Using (2.4), we have: 



(2.7) 



-z + V- — z 



V 



ipA 



^ ^ (t)A (t>D 



(PD , (pA -1 
<PA (Pd 



(2.8) 



M 



( (l)C (1>B _i 

0B 0c 



0s "/"cA \ 
4>c <Pb 



V 



<PB (pC -1 

(PC 4>B 



(2.9) 



Discrete zero curvature condition for these operators has the form (1.9). It gives rise to a system 
of non-hnear equations for fields at the vertices of four neighbouring elementary cells, which is a direct 
corollary of the basic equation (2.1). 

Let us turn to the lattice SG model with continuous time introduced by Izergin and Korepin in |^]. 
Its L-operator on l-th site has the form[| 



zxi +z ^Xi^ 



S ^ifilTTl 



zxi ^ + z ^xi 



(2.10) 



Here ttj, xi ^-re exponentiated canonical variables on the lattice, ipi ^ [l + s{xf + Xi ^)] ^ ^^J^d s is a 
parameter of the model. 

As is known, discretizing time in integrable models consists merely in changing the M-operator and 
leaving the L-operator unchanged, so we are to identify the L-operators ( 2.10| ) and (2.8). For that 
purpose, consider composite fields 



7r(u, v) = [ip{u + 1, v)^l^{u, V + , x(u, v) 



(2.11) 



on the chiral lattice and set tt/ = 7r(/, I), xi = xih ^-t the constant time slice to = 0. Finally, identifying 
s = /ii^(/i^ + v'^)~^ and taking into account the equation of motion (p.l[), we conclude that 



'.{IK) 



{z)^iH Ci{{^iv)'^ z) . 



(2.12) 



Here we use the natural notation 

Ci{z) = Cd.^aAz) , (2.13) 

where Ai — {I, I), Di = {I + 1,1 + 1) (note that in this notation Di — Ai+i). (In what follows we use the 
similar notation for the Af-operator (2.9): A^s,^^, (^) = Aii{z), where Bi — {I + 1,1 — 1).) 



The L-operator c'l'^^\z) has two degeneracy points z^. In terms of the parameters ^,1^ they are 
z^ — {^/ly)^^. Using ( 2.12 ) and (2_^), it is not difficult to represent £|^^''(z^) in the form ( |l.4| ) with 
the r.h.s. expressed through the field ili{u,v). 



3 Hirota's bilinear formalism 



The idea of Hirota's approach is to treat eq. (2.1) as a consequence of 3-term bilinear equations for 
T-functions, which are simpler and in a sense more fundamental. In this section we give a minimal 
extraction from Hirota's method necessary for what follows (for more details see e.g. jl^, |p^). 

^Wc use a slightly modified form of the L-operator originally suggested in They differ by a unitary transformation 
and multiplication by the mat rix m from the left. The latter modification takes into account the fact that we deal with 
the Faddeev-Volkov equation (|2.l|) rather than Hirota's discrete SG equation itself. 
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In the case at hand we need two r-functions: t and f. Making the substitution 

, , , t(u, v) ,„ , 

^P{u,v)^^-^ (3.1) 
r(w, v) 



we immediately see that eq. (2.1) fohows from the couple of Hirota's bilinear equations g 

{v - h)taTd = vtbTc - ^.tbTc , 

{v - iJ.)TAfD = vtbTc - HTbTc . (3.2) 
The equivalent form of these equations, 

{V + pL)TBfc = IJ-TATD + VTaTd , 

{v + ij.)tbtc = IJ-TATD + vtatd , (3.3) 
is equally useful in what follows. At last, we point out the relation 

t{u — 1, v)t{u + 1, u) + t{u — 1, v)t{u + \,v) — 2t(u, w)f (u, v) (3.4) 
valid for the r-functions of the discrete SG model. 



Let us give some additional remarks. Equations (3.2) are to be thought of as being embedded into 
the infinite 2D Toda lattice hierarchy with 2-periodic reduction jl^. This embedding implies that 
parameters /i and y are Miwa's variables corresponding to the two elementary discrete flows 
Miwa's variables play the role of inverse lattice spacings in the chiral directions. Lattice spacing in the 
m-direction is then to be identified with ^^j^. One should note, however, that in this approach the chiral 
coordinate axes are in general not orthogonal to each other. In particular, as it is seen from eqs. ( |3.2|) , 
&t ^ — u one must identify u and v, so that the 2D lattice collapses to a ID one. In this sense eq. (3.4) 
follows from (3.3) ai v ~ ^. 



4 Quantum i?-matrix in the classical lattice SG model 



Our goal is to represent the L and M operators (| 
some vectors in the "quantum" space. 

Consider quantum i?-matrices of the form 

/ a{z) 



, ( 2^ ) as convolutions of quantum i?- matrix with 











\ 



±b{z) c{z) 
c{z) ±h{z) 



V 











(4.1) 



i{z) ) 



where 



a(z) = qz — q 



b{z)=z-z-\ c{z) = q-q-\ (4.2) 

Here q is a "quantum" parameter and z is the spectral parameter. The i?-matrices _R^+^ and i?^^) differ 
by Drinfeld's twist. Both of them satisfy the quantum Yang-Baxter equation. In the Introduction we 
have used the notation _R(+)(z) = R{z). When necessary, we write R{z) = R{z;q). 

The i?-niatrix has 2x2 block structure with respect to the "auxiliary" space. Let i,i' number block 
rows and columns and let number rows and columns inside each block (i.e., in the "quantum" 



space). Then matrix elements of the i?-matrix of the form (4.1) are denoted by R{z)Jj,. Let |a), |/3) 



be two vectors in the quantum space (see ( |1.5| )). Each block of the _R- matrix is an operator in the 
quantum space. Consider its action to |a) and subsequent scalar product with (see the figure in the 
Introduction). The result is a 2x2 matrix in the auxiliary space: 



{f3\R{z)\ 
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Substituting the matrices (4.1), we find: 

/3(i)a(i)a(z) ± /3(2)a(2)5(^) /3(2)q,(i)c(z) 

{/3\R'-^\z)\a) = ' 

/3(i)a(2)c(z) ±/3(i)a(i)6(z) + /3(2)a(2)a(z) 



(4.3) 



We are going to identify r.h.s. of eqs. ( p.8[ ), (2.9) with this matrix. The key step is to express matrix 
elements of the L and M operators thro ugh the t- functions according to (3J) and after that make use 
of Hirota's bilinear equations ( |3.2[ ), (3.3) when necessary. The best result is achieved after the diagonal 
gauge transformation 



CD^Aiz) — > Ca^dIz) = ( ^ J^ ] Cd^a{z) , 

\ TATA ) 



\TCTCJ 



Mb^c{z) - 

Omitting details of the computation, we present the final result. Set 



\P) = 



q = 



Using eqs. (3.2), (|3.3|), we get 



(/3z5|i?W(z)|aA) = 



11 + v 



TATACo^AitJ-z) , 

TcTcMB^-cifJ-z) ■ 



(4.4) 
(4.5) 

(4.6) 

(4.7) 
(4.8) 



Finally, let us pass to the discret e spa ce-time coordinates /, m and consider the slice of the lattice at 
TO = 0. The notations of the type ( 2.13 ) explained in Sect. 2 are again convenient. Making use of eq. (^^), 
we arrive at the formulas 



Miifiz) = 



2^iu {l3i\R^-\z)\ai) 



(4.9) 



(4.10) 



which up to the constant prefactors coincide with the ones announced in the Introduction. Location of 
the vectors |q;;), 1/3;) on the lattice with respect to the m = slice is shown in the figure: 




Explicitly, they are: 



\ai) 



f(/ + l,0 

T{l + l,l) 



The normalization factor in eq. (1.4) is equal to A; = ^Jiv{^ — v) ^t{1, 1)t{1, I) 




(4.11) 
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5 On the continuous time limit 



In this section we show that the r-matrix formula for local Af-operators ( |1.6| ) is a degenerate case of 



eq. (1.8). As it was already mentioned at the end of Sect. 3, plays the role of lattice spacing of the 



discrete time variable m. A naive continuous time limit would then he v ^ that means g — s- 1 in the 
i?-matrix, so the first non-trivial term in its expansion in g — 1 is just the classical r-matrix. This agrees 



with eq. (1.6). However, this limiting procedure would imply linig^i = that is certainly wrong 
in general (see e.g. the figure in the previous section). The correct limiting transition to a continuous 
time coordinate needs some clarification. 

The limit — > changes parameters of the space lattice and the L-operator. That is the reason why 
the naive limit does not work. In the correct limiting procedure, the time lattice spacing must approach 
zero independently of and v. 

To achieve the goal, we introduce v' - an other "copy" of the chiral flow v with the parameter i/', 
so now we have a 3D lattice. Equations ( |3.2| ) are valid in each 2D section of this lattice of the form 
v' — const, u = const oi v — const. For instance, in the latter case we have (cf. (3.2)): 



(h'' — ll)f{u,V,v')T(u+l,V,v' + l) = l''T{u+l,V,v')f{u,V,v' + 1) — j.lT{u+l ,V ,v')t{u,V ,v' + 1) , 

— ^)t{u,v,v')t{u+1,v,v' + l) = h''T{u + l,v,v')T{u,v,v' + 1) — fiT{u+l,v,v')f{u,v,v' + 1) . (5.1) 
Now we can tend v' ^ leaving unchanged. Set 

q' = Ii = l + £ + 0{e^), e^O. (5.2) 

The small parameter e has the meaning of lattice spacing in the auxiliary direction ni' ~ ^{u — v'). 

Discrete A/-operators are defined up to multiplication by a scalar function of z independent of 
dynamical variables. To pursue the continuous time limit, it is convenient to normalize the Af-operator 
in such a way that Aii(z) = J at e = 0. Then the next term (of order e) yields a local Af-operator of 
continuous time fiow. To find such an Af-operator at the l-th site of the ID lattice (diagonally embedded 
into the u, w-lattice as {I, I)), we should expand the discrete Af-operator A^^j'^c,' {z) which generates the 
shift {I — 1, 1, 1) — > {1, 1, 0) on the 3D lattice with coordinates (u, v, v'). This is illustrated by the figure 




which shows the u,w'-section of the 3D lattice spanned by u,v,v'. Coordinates of the parallelogram 
vertices are: A[ = {I - 1, 1, 0), = Ai = {1, 1, 0), C/ = {I - 1, 1, 1), = {1, 1, 1). The point C/ tends to 
the point B'l = Ai as i^' — > ^, so the parallelogram collapses to the u-axis. We have: 



where 



MiipLz) 



MB[^c[{z)^I + eMi{z) + 0{e''), 

t{1-1,1)t{1 + 1,1) 



/I _,^ t{1-1,1)t{1 + 1,1) 

2^ ' r{l,m,l) 

f{l - 1,1)t{1 + 1,0 

Til,l)T{l,l) 



(5.3) 



t{1,1)t{1,1) 

1 _i^ f(/-l,OT(/ + l,0 

2^" + " ^ r{l,l)f{l,l) 



(5.4) 
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By definition, the classical r-matrix is 

R^+\z-q') - (z - z-^)I(g,I 



r{z) = lim 



e^o e{z — z^i) 

[(-^ + ■^"^)-^ / + 2cti ® (Ti + 2(72 (72 + (Z + 2:"^)(73 (g) (73] . (5.5) 



Comparing with (5.4), we obtain the r-matrix representation 

MA^^) = (5.6) 



{Hi 


r(z) 


OLl-\) 


(A 




-l) 



of the Af -operator (see (pT^)) with the trigonometric classical r-matrix (5.5) 



6 Conclusion 



The main result of this work is the i?-matrix representation (4^), ( 4.10| ) of the local L-M pair for the 



classical lattice SG model with discrete space-time coordinates. In our opinion, the very fact that the 
typical quantum i?-matrix naturally arises in a purely classical problem is important and interesting 
by itself. As a by-product, we have shown that components of the vectors |a), (representing the 



L-operator at the degeneracy point, see (1.4)) are r-functions. 

Let us recall that the quantum Yang-Baxter equation already appeared in connection with purely 
classical problems, though in a different context |2^. However, the class of solutions relevant to 
classical problems is most likely very far from i?-matrices of known quantum integrable models. In our 
construction, the role of the quantum Yang-Baxter equation remains obscure; instead, the most popular 
4x4 trigonometric quantum i?-matrix is shown to take part in the zero curvature representation of the 
classical discrete SG model. We believe that a conceptual explanation of this phenomenon nevertheless 
relies on the Yang-Baxter equation. 

We should stress that the " quantum deformation parameter" q of the trigonometric i?-matrix in our 
context seems to have nothing to do with any kind of quantization. Indeed, it is related to the mass 
parameter and the lattice spacing of the classical model. This fact suggests to ask for a non-standard 
hidden i?-matrix structure of the model, which might survive and show up on the quantum level, too. 



Another interesting question is to find a discrete time analogue of the non-local generating function (1.3) 
of M-operators. 

At last, let us point out that the i?-matrix representation can be extended to il/-operators of a 
more general model - partially anisotropic classical Heisenberg spin chain in discrete time. This system 
includes the discrete KdV and some other equations as particular cases. The results will be reported 
elsewhere. 
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